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This  report  presents  the  results  of  an  experimental  investiga¬ 
tion  of  the  dynamic  stall  of  two  rotor  blade  airfoil  sections 
(NACA  0006  and  Vertol  13006-.7)  chosen  for  high  Mach  number 
operation.  Forces  and  moments  in  two-dimensional  flow  were 
determined  by  measuring  differential  pressures  in  pitching 
oscillation.  Mach  number,  Reynolds  number,  and  mean  angle-of- 
attack  ranges  tested  corresponded  to  the  conditions  prevailing 
on  the  retreating  blades  of  full-scale  helicopters.  Frequen¬ 
cies  of  oscillation  ranged  up  to  the  level  corresponding  to 
typic-’’'  first  torsional  vibration  modes  of  rotor  blades. 

Pitch  oscillation  was  Tound  to  increase  the  angle  of  attack  at 
which  stall  occurred,  with  a  corresponding  increase  in  maxi¬ 
mum  normal  force.  Hysteresis  effects  on  the  pitching  moment 
response  produced  regions  of  negative  aerodynamic  damping  in 
oscillations  at  mean  angles  of  attack  near  the  steady-flow 
stall  values  for  all  Mach  numbers  from  0.2  to  0.6. 

The  reduced  frequency  of  oscillation  was  found  to  be  the  domi¬ 
nant  influence  on  dynamic  stall  behavior,  although  some  com¬ 
pressibility  effects  were  noted  at  Mach  numbers  up  to  0.4,  and 
strong  qualitative  changes  in  the  stall  behavior  were  evident 
beyond  that  value. 

Leading-edge  camber  and  increased  nose  radius  were  found  to 
increase  the  angle  of  attack  and  normal  force  at  stall  at  all 
frequencies  tested,  as  well  as  at  fixed  angle  of  attack.  Pat¬ 
terns  evident  in  the  force  and  moment  data  indicate  that  these 
design  features  produce  an  effect  much  like  a  reduction  of  the 
Mach  number. 
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FOREWORD 


The  results  of  a  program  of  testing  thin  airfoils  oscillating 
near  stall  are  summarized  in  this  report.  The  project  was  per¬ 
formed  for  the  United  States  Army  Aviation  Materiel  Laborator¬ 
ies  under  Contract  DAAJ02-67-C-0095  and  under  the  technical 
cognizance  of  Clifton  G.  Wrestler  and  Patrick  Cancro  of  the 
Aeromechanics  Division  of  USAAVLABS. 

The  report  consists  of  two  volumes: 

Volume  I,  Summary  and  Evaluation  of  Results 

Volume  II «  Data  Report 

The  tests  were  conducted  at  the  Commercial  Airplane  Division 
of  The  Boeing  Company  in  the  supersonic  wind  tunnel.  The 
assistance  and  cooperation  of  the  Model  Design,  Instrumenta¬ 
tion,  and  Supersonic  Wind  Tunnel  Testing  Groups  are  gratefully 
acknowledged.  The  assistance  of  Messrs.  W.  McLean,  K.  Simmonds, 
and  L.  Dadone  of  the  Vertol  Division  in  running  the  tests  is 
also  greatly  appreciated. 
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INTRODUCTION 


The  speed  and  lifting  capability  of  contemporary  helicopters 
are  limited  by  retreating  blade  stall.  As  aircraft  speed  or 
weight  is  increased  beyond  the  stall  limit,  rapid  increases 
in  blade  loads,  control  loads,  or  rotor  horsepower  occur. 
Attempts  to  predict  these  effects  by  using  steady-flow  airfoil 
characteristics*  in  strip  theory  analyses  have  been  unsuccess¬ 
ful;  the  onset  of  stall  is  predicted  too  early. 

There  are  two  recognized  aerodynamic  effects,  apart  from  cen¬ 
trifugal  effects  on  the  boundary  layer,  that  can  cause  the 
differences  between  test  and  theory:  yawed  flow  and  unsteady 
aeiodyncunics .  Neither  effect  has  been  included  in  current 
rotor  analyses.  In  discussing  the  effects  of  yawed  flow  on  the 
rotor  blade,  Harris^  shows  that  part  of  the  discrepancy  between 
test  and  theory  can  be  eliminated  by  including  yaw  effects  in 
the  airfoil  tables.  The  second  effect,  unsteady  aerodynamics, 
has  been  shown  by  various  investigators  (see  Bibliography  in 
Reference  2)  to  produce  a  stall  delay  which  increases  the  max¬ 
imum  lift  coefficient  above  that  for  static  stall.  In  addition 
to  this  stall  delay,  negative  aerodynamic  deunping  can  occur  for 
pitch  oscillations  through  stall.  This  negative  damping  is 
caused  by  time  lag  effects  in  the  blade  pitching  moment  versus 
angle-of-attack  relation  and  can  cause  large  torsional  blade 
deflections  and  large  control  loads.  The  resulting  blade  and 
control  stresses  caused  by  stall  can  limit  the  flight  speed  of 
the  helicopter;  therefore,  it  is  important  that  the  effects  of 
dynamic  stall  be  fully  understood.  Further,  these  dynamic 
effects  must  be  incorporated  in  rotor  blade  analyses  in  order 
to  predict  the  performance  and  flight  limitations  of  high- 
performance  helicopter  rotors  noar  stall. 

A  comprehensive  literature  search  revealed  little  useful  oscil¬ 
latory  data  appropriate  to  helicopter  blade  applications  and 
no  theories  which  adequately  predict  dynamic  stall  effects. 

To  establish  a  comprehensive  framework  of  experimetnal  data  for 
immediate  use  in  rotor  analyses,  and  to  provide  sufficient 
data  to  stimulate  the  formulation  of  a  dynamic  stall  theory, 
a  program  of  two-dimensional  wind  tunnel  tests  was  conducted 
in  1966-1967  and  is  reported  in  Reference  2.  The  two  airfoils 
tested  in  that  program  were  the  NACA  0012  (modified)  used  on 
the  CH-47A,  and  the  Vertol  23010-1.58  currently  used  on  the 
CH-47B  and  C  models,  xhose  tests  substantiated  the  work  of 
previous  investigations  and  carried  the  results  into  many  areas 
where  no  data  were  formerly  available. 


^Tabulated  aerodynamic  coefficients  reflecting  both  Mach  num¬ 
ber  and  flow  separation  effects. 
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Many  current  high-performance  rotor  blade  designs  use  thin 
airfoil  sections  near  the  tip  to  postpone  the  onset  of  power 
losses  due  to  compressibility. 

No  dynamic  stall  data  for  such  profiles  were  available,  and 
the  need  was  clear.  The  test  program  reported  here  was  con¬ 
ducted  to  fill  that  need. 
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TEST  FACILITIES,  MODELS.  AND  DATA  SYSTEM 

The  section  presents  a  brief  description  of  the  basic  test 
facility,  apparatus,  and  data  recording  and  reduction  system. 
These  are  described  in  more  detail  in  Reference  2.  For  the 
present  tests,  a  number  of  improvements  (relating  principally 
to  model  construction  techniques  and  to  ^e  data  reduction 
system)  have  been  made  and  are  discussed  below. 

TEST  FACILITIES 

The  subsonic,  variable-density,  two-dimensional,  l-by-3-foot 
test  section  of  the  Boeing  4-by-4-foot  supersonic  wind  tunnel 
was  used  to  test  8- inch-chord  model  airfoils.  The  test  section's 
total  pressure  and  airfoil  drive  frequency  were  increased  to 
simulate  Reynolds  numbers  and  reduced  frequenices  correspond¬ 
ing  to  a  2-foot-chord  full-scale  rotor  blade.  The  tunnel  Mach 
number  capability  with  the  model  installed  was  limited  to  0.75 
because  of  blockage  effects. 

The  pitch  oscillating  mechanism  is  shown  in  Figure  1.  A 
hydraulic  motor  supplies  the  power  to  rotate  a  flywheel  with 
an  eccentrically  mounted  cam.  The  cam  slides  in  a  slotted 
crank  which  is  rigidly  attached  to  the  airfoil  model  on  the 
pitch  axis.  Four  fixed  nominal  eunplitudes  of  oscillation, 
from  2.5  to  10  degrees,  are  provided  by  different  eccentric 
locations  for  the  c^ufn. 

AIRFOIL  MODELS  AND  INSTRUMENTATION 

The  first  test  airfoil  was  the  Vertol  13006-.7.  This  nose- 
ceunbered  airfoil,  of  6  percent  thickness  ratio,  has  an  in¬ 
creased  leading-edge  radius  compared  with  the  stemdard  NACA 
13006  section.  This  modification  improves  the  maximum  lift 
characteristics  without  raising  the  drag  level.  The  aerodynam¬ 
ic  design  considerations  of  this  airfoil  section  are  reported 
by  Davenport  and  Front^.  To  provide  a  broader  data  base  for 
thin  airfoil  dyneunic  characteristics,  a  thin  symmetrical  air¬ 
foil  section  (the  NACA  0006)  was  also  tested. 

Airfoil  contours  and  coordinates  for  each  model  are  shown  in 
Figure  2;  a  statement  of  the  basic  model  parameters  and  instru¬ 
mentation  Installation  is  given  in  the  table  on  page  6.  Beth 
models  were  machined  from  solid  steel  in  order  to  obtain  ade¬ 
quate  structural  stiffness. 

The  miniature  differential  pressure  transducers  used  for  this 
test  were  encapsulated  in  cylindrical  plugs  before  installation 
in  the  airfoil,  as  shown  in  Figure  3. 


Pitch  Oscillating  Mechanism. 
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Figure  2.  Test  Airfoil  Coordinates. 
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AIRFOIL  MODEL  AND  INSTRUMENTATION  DATA 


I  tern 

Description 

Airfoil  Section 

NACA  006  and  Vertol  13006- . 

7 

Model  Span 

12  inches  nominal 

Model  Scale 

1/3,  based  on  2- foot-chord 
rotor  blade 

Model  Chord 

8  inches 

Thickness  Ratio 

6  percent 

Construction 

Machined  from  maraging  steel 

Transducer  Type 

Scientific  Advances  SD-M6 

Number  Installed 

14  per  model 

Location  in  Percent 

Chord 

1,  2.5,  5,  10,  15,  20,  25, 
40,  50,  70,  80,  87.5,  95 

30, 

Pressure  Range 

+50  psi  located  from  1  to 
percent  chord 

10 

+  25  psi  located  from  15  to 
percent  chord 

95 

Minimum  Natural  Frequency 
of  Installed  Transducers 

4,000  Hertz 

Damping  Ratio  of 

Installed  Transducers 

0.25 
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DATA  RECORDING  AND  REDUCTION  SYSTEM 


The  complete  data  recording  and  reduction  system  is  shown 
schematically  in  Figure  4.  The  data  recording,  on-line  oscil¬ 
lograph,  and  X-Y  plotter  presentation  were  controlled  automat¬ 
ically  during  test  runs  through  the  tunnel  operating  console. 
After  analysis  and  editing  of  the  oscillograph  records,  all 
further  data  reduction  was  done  off-line  at  the  Vertol  Divi¬ 
sion. 

Data  Recording 

All  test  data,  including  both  airfoil  pressures  and  tunnel 
operating  condition  information,  were  recorded  simultaneously 
by  a  pair  of  Sangamo  Model  3500  FM  wideband,  14-channel  tape 
recorders.  A  time  code  signal  and  a  flywheel  1/rev  pulse  were 
duplicated  on  both  tapes  to  provide  data  synchronization.  In 
addition,  an  Electronic  Associates  Model  TR48  analog  computer 
was  used  to  integrate  the  differential  pressures  to  provide 
preliminary  on-lir values  of  Cn  and  Cn,  which  were  recorded 
for  reference. 

All  test  parameters,  including  Cn  and  Cm/  were  recorded  simul¬ 
taneously  on  oscillograph  stripouts  for  direct  data  monitoring. 
A  time  integral  value  representing  the  average  aerodynamic 
cycle  damping  given  by  ^  CM<da  was  also  included.  In  addition, 
an  X-Y  pen  plotter  was  used  on-line  to  record  scaled  aerodynam¬ 
ic  damping  versus  Oo»  or  Cn  versus  a#  during  steady-flow  tests. 
This  information  formed  the  basis  for  planning  the  detailed 
execution  of  a  major  part  of  the  test  program. 

Data  Reduction 


Six  hundred  and  seventy-three  different  test  conditions  were 
recorded  and  digitally  processed,  representing  a  total  of  more 
than  8  million  digital  readings.  Each  oscillatory  test  point 
represented  an  average  of  at  least  10  individual  consecutive 
data  cycles  (20  cycles  at  drive  frequencies  of  48  Hertz  and 
higher)  ,  with  approximately  50  points  read  for  each  cycle  and 
parameter. 

To  process  this  large  volume  of  data,  the  off-line  data  reduc¬ 
tion  system  described  in  Reference  2  was  used  to  convert  the 
analog  signals  to  digital  form  and  to  provide  all  final  data 
output.  Two  additional  corrections  were  made  to  the  data  taken 
from  these  tests  to  account  for  the  significant  nonlinear  sen¬ 
sitivity  characteristics  of  the  transducer  units  and  their 
response  to  the  pitch  accelerations  of  the  model. 
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EXPERIMENTAL  RESULTS 


PARAMETRIC  TRENDS  FOR  DYNAMIC  Cn  AND  Cm  BEHAVIOR 


Forced  pitch  oscillation  test  data  were  taken  for  the  symmetri¬ 
cal  NACA  0006  and  the  cambered  Vertol  13006-.7  airfoil  sections. 
Nominal  values  of  mean  angle  of  attack  ranged  from  0  to  25 
degrees,  eunplitude  from  2.5  to  10  degrees,  Mach  number  from  0.2 
to  0.75,  and  reduced  frequency  from  0.04  to  0.66.  The  experi¬ 
mental  results  are  usually  discussed  and  identified  in  terms  of 
selected  nominal  values  of  these  parameters.  Exact  values  of 
the  parameters  are  used  in  the  discussion  only  when  required  to 
identify  specific  items.  However,  as  with  any  test  program,  the 
exact  values  of  the  parameters  varied  slightly  from  test  point 
to  test  point.  Therefore,  exact  values  of  M,  a^,  Ao,  f,  and  k 
are  shown  on  each  plot  where  applicable.  The  variations  are  not 
sufficiently  large  to  require  cross-plotting  of  the  data  to 
determine  trends  and  levels  at  the  nominal  conditions . 

Figures  5  through  10  present  the  effects  of  flow  and  oscillation 
par2uneters  on  the  Cn  and  Cm  versus  a  behavior  for  the  Vertol 
13006-.7  only,  since  the  NACA  0006  data  are  similar  except  where 
specifically  noted.  The  aerodynamic  characteristics  for  both 
airfoils  are  also  similar  to  the  results  obtained  for  the 
thicker  airfoils  of  Reference  2. 

Effect  of  Mean  Angle  of  Attack 

Figure  5  presents  Cn  and  Cm  versus  a  traces  for  a  5 -degree  amp¬ 
litude,  low-frequency  oscillation  of  the  Vertol  13006-.7  air¬ 
foil  at  M  a  0.2  and  k  «  0.11.  A  sequence  of  Oo's,  from  5  to 
17.5  degrees,  shows  the  development  of  Cn  and  Cm  from  fully 
attached  flow  conditions  to  those  with  fully  stalled  flow. 
Steady-flow  data  are  included  on  the  plots  for  reference. 

The  first  pair  of  traces,  at  =  4.93  degrees,  shows  no  sign 

of  stall.  The  Cn  and  Cm  curves  are  of  the  characteristic 
elliptical  form  predicted  by  potential- flow  theory.  (Both  the 
Cn  and  the  C^  curves  are  displaced  from  the  steady-a  lines. 

This  is  not  considered  to  be  aerodynamics lly  significant,  and 
is  attributed  to  baseline  pressure  measurement  error  due  to 
the  low  dyneunic  pressure  at  this  Mach  number.)  The  next  pair, 
at  Oo  7.48  degrees,  shows  the  first  signs  of  stall  in  the 
figure-eight  shapes  of  both  Cn  and  Cm>  Note  the  delayed  stall 
and  the  increases  in  maximum  normal  force  and  negative  moment 
peak  compared  to  the  static  values.  The  stall  angle  goes  up 
as  Oe  is  increased  further,  but  it  occurs  increasingly  earlier 
in  the  cycle.  This  is  associated  with  further  gains  in  the 
maximum  Cn>  which  reaches  a  peak  for  Qq  =  14.82  degrees.  When 
Oo  is  increased  to  17.31  degrees,  the  minimum  angle  of  attack 
during  oscillation  is  now  above  the  static  stall  value  and  is 
too  high  for  there  to  be  a  complete  recovery  to  attached  flow. 
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characteristics.  Consequently,  there  are  reductions  in  meuci- 
mum  Cn  and  in  the  amplitudes  of  the  Cn  and  Cm  excursions, 
although  they  remain  approximately  centered  on  the  steady  data. 

Effect  of  Pitching  Frequency 

Figure  6  shows  a  sequential  set  of  Cm  emd  Cm  traces  at  increas¬ 
ing  frequency  for  the  Vertol  13006-.7  airfoil  at  M  =  0.2,  La  ■ 

5  degrees,  and  Uq  *  10  degrees.  The  most  visible  trend  with 
frequency  is  the  increased  postponement  of  stall  effects, 
evidenced  by  the  overshoot  behavior  of  both  Cn  and  Cm  in  com¬ 
parison  to  the  superposed  steady-flow  characteristics. 

The  Cm  break  at  stall  occurs  initially  during  the  increaslng-a 
part  of  the  cycle.  The  stall  effect  is  delayed  as  frequency 
is  increased  until,  at  k  «  0.351,  it  occurs  at  the  point  of 
maximum  a  in  the  cycle.  Further  increases  in  frequency  result 
in  a  progressive  postponement  into  the  decreasing-a  portion  of 
the  cycle.  Potential-flow  Cm  ellipses  at  the  lowest  and  high¬ 
est  frequencies  are  given  for  comparison. 

The  area  enclosed  by  the  Cm  trace  is  proportional  to  the  work 
per  cycle  (or  "cycle  damping")  and  has  importeuit  physical 
significance.  The  net  work  done  by  the  airfoil  on  the  surround¬ 
ing  air  is  proportional  to  the  integral 


W  »  Cm  da.  (1) 


The  damping  is  positive  for  a  counterclockwise  circuit.  Areas 
enclosed  by  a  clockwise  circuit  are  damped  negatively;  i.e., 
the  airfoil  system  extracts  energy  from  the  airstream.  This 
can  lead  to  an  increase  in  the  amplitude  of  oscillation  with 
time  for  an  elastic  system,  and  it  is  precisely  the  condition 
for  flutter. 

The  Cm  traces  demonstrate  the  effect  of  stall  delay  on  the 
aerodynamic  damping.  As  the  a  for  Cm  stall  increases  to  the 
maximum  value,  the  net  area  contributing  to  negative  damping 
also  grows  to  a  maximum.  At  higher  frequencies,  the  Cm  break 
occurs  after  maximum  a.  This  creates  a  counterclockwise  end 
loop,  whose  growth  ultimately  leads  to  net  positive  damping. 

Up  to  the  frequency  where  the  Cm  break  occurs  at  the  maximum  a, 
^MAX  occurs  after  the  Cm  break  and  coincides  approximately 
with  the  negative  moment  peak.  At  higher  frequencies,  the 
Cnmax  position  reverts  to  a  point  where  a  is  still  increasing. 
This  is  consistent  with  potential  flow  theory,  as  shown  for 
the  case  at  the  highest  frequency. 
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Note  that  Cn  tends  briefly  to  move  toward  a  more  positive  value 
at  the  moment  when  Cm  breaks  in  the  negative  direction.  This 
behavior,  which  will  be  discussed  later  in  terms  of  pressure 
distribution  data,  indicates  that  a  vortex  is  shed  from  the 
leading  edge  at  the  moment  of  flow  separation. 

Effect  of  Pitch  Amplitude,  Maximum  a  Fixed 

The  Cn  and  Cm  traces  for  three  different  amplitudes  of  oscilla¬ 
tion  at  M  =  0.2  and  k  =  0.11  compare  the  stall  behavior  at 
constant  maximum  a  in  Figure  7.  Since  the  nominal  frequency 
of  oscillation  is  the  same  for  each  case,  the  differences  in 
the  Cn  and  Cm  characteristics  are  the  result  of  pitch  rate 
effects,  da/dt  being  proportional  to  Aa. 

Cnmax  increases  noticecibly  with  increasing  amplitude.  In  addi¬ 
tion,  the  angles  of  attack  for  Cm  stall  break  and  for  CnmaX 
higher  at  increasing  amplitude  and  do  not  correspond  to  a  fixed 
a  increment  above  the  static  stall  value.  These  effects  are 
similar  to  those  induced  by  frequency  increases,  as  discussed 
earlier. 

Effect  of  Mach  Number 

To  illustrate  the  effects  of  Mach  number  on  stalling  character¬ 
istics,  Figure  8  presents  three  sets  of  Cn  and  Cm  traces  at 
different  parametric  conditions  for  Mach  numbers  of  0.2  and 

0.4. 

The  firct  set  of  traces,  at  f  =  12  Hertz  (scale  1/rev  frequen¬ 
cy)  t  Oo  =  10  degrees,  and  Aa  =  5  degrees,  shows  reductions  in 
the  Cnmax  stall  with  an  increase  in  Mach  number.  In 

addition,  the  Cn  and  Cm  loops  are  smaller.  These  changes, 
however,  must  be  ascribed  to  the  effect  of  the  reduction  of 
reduced  frequency  (k)  by  inference  from  the  second  and  third 
sets  of  curves. 

The  second  set  of  curves  illustrates  the  comparative  behavior 
with  Mach  number  at  nearly  constant  reduced  frequency.  In  this 
case,  close  similarities  in  the  maximum  excursions  of  the  Cn 
and  Cm  behavior  are  evident.  However,  the  pitching  moment 
traces  do  indicate  an  earlier  inception  of  stall  at  the  higher 
Mach  number  condition,  which  substantially  reduces  the  net  area 
of  the  Cm  loop.  This  is  coupled  with  only  small  changes  in 
the  Cn  trace  close  to  the  position  of  Cnmax* 

A  comparison  of  Mach  number  effects  at  the  same  k  but  with  a 
better  match  of  the  stalling  behavior  is  presented  in  the  third 
set  of  traces.  Here,  Oo  is  reduced  by  2.5  degrees  for  the 
M  =  0.4  case.  The  traces  show  a  remarkable  similarity  in  shape, 
which  emphasizes  the  sensitivity  of  the  stalling  to  Mach  number, 
allowing  for  the  expected  displacement  due  to  the  difference 
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in  ao.  (These  strong  similarities  do  not  persist  much  beyond 
M  =  0.4,  because  of  the  influence  of  transonic  flow  on  the 
stall  process,  as  will  be  made  clear  later.  For  Cn  and  Cm 
traces  at  M  =  0.6,  the  reader  is  referred  to  Volume  II,  Figure 
1.) 

Effect  of  Leading-Edge  Design 

Figure  9  presents  a  comparison  of  the  dynamic  Cn  and  Cm  behavior 
between  the  Vertol  13006-.7  and  NACA  0006  airfoils  at  M  =  0.2 
and  f  =  12  Hertz  (equivalent  1/rev  frequency) .  Note  that  the 
a  for  Cm  stall  break  is  about  3  degrees  higher  for  the  Vertol 
13006-.7  airfoil,  which  corresponds  closely  to  the  difference 
in  the  steady-flow  a  for  the  inception  of  laminar  separation. 
(See  Appendix  I  for  discussion.) 

To  illustrate  the  effect  of  camber  and  leading-edge  radius  at 
higher  frequency  and  Mach  number.  Figure  10  presents  Cn  and 
Cm  loops  at  M  =  0.4  and  f  =  48  Hertz.  Because  a  comparison  of 
the  two  airfoils  at  the  same  oiq  would  indicate  significant 
dissimilarities  in  the  traces,  the  ao  has  been  reduced  by  2.2 
degrees  for  the  NACA  0006.  The  similarities  between  the  stall¬ 
ing  patterns  for  the  airfoils  are  striking.  The  difference  in 
ao  is  consistent  with  steady-flow  stalling  behavior  presented 
in  Appendix  I,  from  which  it  is  inferred  that  airfoil  geometry 
effects  observed  in  steady  flow  persist  under  oscillatory  con¬ 
ditions  . 

Note  also  that  the  traces  in  Figure  10  show  a  remarkable  simi¬ 
larity  to  those  illustrating  Mach  number  effects  presented  in 
Figure  8.  The  most  important  difference  among  the  traces  is 
the  value  of  ao- 

SUMMARY  OF  CNj^j^  AND  DAMPING  DATA 


Maximum  Dynamic  Normal  Force 

Figures  11  through  13  show  the  maximum  lift  capabilities  of  the 
NACA  0006  and  Vertol  13006-.7  airfoils.  For  a  realistic  com¬ 
parison  of  the  lift  potential  of  these  two  airfoils  with 
respect  to  the  flight  operating  environment  of  a  rotor  blade, 
it  is  necessary  to  study  behavior  in  1/rev  motion.  Figure  11 
shows  the  dynamic  maximum  Cn  and  negative  Cm  attained  at  Aa  = 

5  degrees  and  f  =  12  Hertz.  The  cambered  airfoi]  shows  sub¬ 
stantially  better  lift  for  all  Mach  numbers  tested  (0.2  to  0.6). 
The  negative  pitching  moment  maxima  are  roughly  equal.  The 
increase  in  dynamic  maxima  with  reduction  of  M  is  due  to  the 
increased  reduced  frequency. 

The  effects  of  the  full  range  of  tost  frequency  and  Mach  number 
on  the  maximum  Cn  values  attained  during  pitch  oscillation  at 
Aa  =  5  degrees  are  shown  for  the  Vertol  13006-.7  airfoil  in 
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Figure  12  and  for  the  NACA  0006  airfoil  in  Figure  13.  The 
reduced  frequency  shows  a  dominant  influence,  while  Mach  number 
effects  are  present,  but  weaker. 

The  effect  of  Mach  number  on  dynamic  stall  is  to  decrease  the 
m^ucimum  attainable  Cm  by  about  the  saune  order  of  magnitude  as 
under  steady-flow  conditions,  except  atM=0.2.  AtM=0.2 

it  was  apparent  that  Oo's*  higher  than  those  tested  here  would 
produce  higher  maximum  Cn's  for  values  of  k  above  0.44;  there¬ 
fore,  these  lines  represent  a  limit  due  only  to  the  scope  of 
testing.  The  data  trend  differences  which  are  apparent  at 
M  »  0.2  are  attributed  to  the  virtual  adosence  of  compressibili¬ 
ty  effects;  here,  the  cambered  airfoil  shows  significantly 
better  maiximum  lift  at  k  =  0.11. 

Aerodynamic  Pitch  Damping 

Principal  pitch  damping  characteristics  for  the  Vertol  13006-.7 
airfoil  are  shown  in  Figures  14  through  16.  Pitch  damping  is 
expressed  as  a  ratio  to  the  theoretical  potential-flow  damping, 
the  derivation  of  which  is  based  on  unsteady  aerodynamic 
theory^  and  has  been  documented  in  Reference  2. 

Figure  14  presents  the  variation  of  pitch  damping  with  fre¬ 
quency  at  M  =  0.2  and  Aa  ^  5  degrees.  Correspondence  with 
potential-flow  damping  is  indicated  by  a  pitch  damping  of 
unity,  and  fair  agreement  is  obtained  until  destabilizing 
effects  due  to  stall  appear.  As  the  frequency  is  increased, 
this  destabilization  is  postponed  to  higher  values  of  Oo, 
as  pointed  out  in  previous  discussion  of  the  parametric  trends. 

Consideration  of  practical  rotor  blade  operation  indicates  that 
the  aerodynamic  damping  characteristics  are  most  important  at 
a  frequency  corresponding  to  the  blade's  first  torsional  mode 
and  at  relatively  small  amplitude.  Figure  15  gives  the  varia¬ 
tion  of  pitch  damping  with  Mach  number  for  f  =  48  Hertz 
(corresponding  to  a  typical  blade  torsional  natural  frequency) 
and  Aa  =  2.5  degrees.  Both  the  amount  of  negative  damping 
and  the  Oo  at  whid'.  it  appears  are  reduced  for  increasing  Mach 
number  until  compressibility  effects  eliminace  negative  damping 
altogether,  as  at  M  =  0.6  and  above. 

Because  of  the  nonlinear  nature  of  stall,  the  damping  can  be 
expected  to  vary  with  amplitude.  This  is  clearly  shown  in 
Figure  16  for  data  at  2.5  and  5-degree  amplitudes  and  f  =  48 
Hertz.  Destabilizing  effects  begin  at  a  higher  Oo  for  Aa  = 

2.5  degrees,  but  the  magnitude  of  the  negative  damping  is 
larger. 

The  mean  angle  of  attack  at  which  the  damping  goes  through  zero 
defines  a  "stability  boundary”,  beyond  which  pitch  oscillations 
may  occur.  Figure  17  presents  stability  boundaries  on  the 
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M  -  Oo  plane  for  the  Vertol  13006-.7  airfoil  to  show  the  effect 
of  2.5  und  5-degree  amplitudes  at  two  frequencies,  correspond¬ 
ing  to  the  rotor  blade  4/rev  and  6/rev  motions.  Mach  number 
effects  limit  areas  of  negative  damping  to  M  <  0.6  at  the  lower 
frequency.  An  increase  in  frequency  leads  to  an  enlargement  of 
the  region  of  instability  and  delays  any  beneficial  compres¬ 
sibility  effects  until  Mach  numbers  greater  than  those  tested. 
Note  that  the  lower  stability  boundary  is  not  strongly  sensi¬ 
tive  to  eunplitude.  Inst2d}ility  generally  appears  at  a  lower 
Qo  for  the  larger  oscillation. 

A  comparison  of  the  stability  limits  for  both  airfoils  at  4/rev 
equivalent  frequency  is  presented  in  Figure  18.  The  NACA  0006 
section  generally  has  stability  boundaries  at  reduced  Oo's  tor 
oscillatory  amplitudes  of  both  2.5  and  5  degrees.  This  is  con¬ 
sistent  with  the  effects  of  leading-edge  design  discussed  pre¬ 
viously.  In  addition,  compressibility  effects  appear  sooner; 
this  is  a  result  of  the  sharper  leading  edge  and  is  substan¬ 
tiated  by  the  steady-flow  results  of  Appendix  I. 

DISCUSSION  OF  STALL  AND  REATTACHMENT  PROCESSES 

Considerable  insight  into  the  flow  processes  governing  dynamic 
Cn,  hysteresis,  and  damping  effects  may  be  derived  from 
study  of  details  of  the  aerodynamic  loading.  To  provide  a 
background  for  comparison  with  the  dynamic  results,  it  is 
essential  first  to  examine  the  variation  of  chordwise  loadings 
experienced  under  steady-flow  conditions. 

NOTE:  The  chordwise  pressure  distribution  plots  presented  in 
Figures  21  through  27  have  been  normalized  by  dividing  ACp  by 
Cn  for  clarity. 

Chordwise  Loadings  Under  Static  Conditions 

Figure  19  indicates  the  progressive  variations  in  the  static 
chordwise  loading  distributions  which  characterize  stall  devel¬ 
opment  on  the  Vertol  13006-.7  airfoil  within  the  range  of  test 
conditions  investigated.  At  M  =  0.2,  the  loading  is  initially 
of  the  classical  potential-flow  type  with  a  high  suction  pres¬ 
sure  peak  at  the  leading  edge,  as  shown  for  a  =  10.2  degrees. 
The  second  loading  distribution,  given  at  a  -  10.6  degrees,  « 
corresponds  to  an  angle  of  attack  just  above  a  small,  sharp 
discontinuity  in  the  Cn  -  a  curve  (see  Figure  36  in  Appendix 
I)  . 

At  a  =  10.6  degrees  the  leading-edge  pressure  peak  shows  a 
significant  reduction,  indicating  that  a  laminar  separation 
bubble  has  formed.  This  phenomenon  is  characteristic  of  "thin 
airfoil  stall"  and  is  discussed  in  detail  by  McCullough  and 
Gault^.  The  sharp  leading  edge  of  a  thin  airfoil  results  in  a 
strong  adverse  pressure  gradient  (just  after  the  peak)  too 


15 


close  to  the  nose  for  transition  from  laminar  to  turbulent 
flow  to  have  occurred  in  the  boundary  layer.  Laminar  boundary 
layers  are  less  resistant  to  separation  than  turbulent  ones, 
and  separation  occurs.  However,  the  very  sharpness  of  the 
pressure  peak  implies  that  the  region  of  strong  adverse  pres¬ 
sure  gradient  must  be  small.  Transition  to  turbulent  separated 
flow  and  reattachment  of  the  boundary  layer  then  result.  The 
gradual  enlargement  of  this  separation  bubble  with  increasing 
a  prevents  further  discontinuities  in  Cn  and  Cm  as  a  increases. 
The  loading  distribution  on  the  fully  stalled  airfoil  (shown 
at  a  =  16.6  degrees)  is  almost  flat,  with  only  a  small  peak  at 
the  leading  edge. 

Next,  consider  the  loadings  at  M  =  0.6,  where  a  very  different 
behavior  is  evident.  Loadings  are  shown  for  angles  of  attack 
at  a  =  8.4  degrees  and  a  =  9.1  degrees,  which  are  just  below 
and  above  the  point  where  deviation  from  the  linear  lift  slope 
occurs  (see  Figure  36  in  Appendix  I) . 

At  a  Mach  number  of  0.6,  the  Cp  corresponding  to  local  sonic 
velocity  is  -1.3.  Assuming  the  most  conservative  pressure 
loading  on  the  lower  surface  (Cp  =  +1.11) ,  the  peak  Mach  number 
on  the  upper  surface  must  be  almost  1.5.  This  implies  that  the 
flow  around  the  nose  is  a  supersonic  expansion,  which  eliminates 
the  strong  adverse  pressure  gradient  noted  for  the  M  =  0.2  con¬ 
dition.  Recovery  to  subsonic  flow  must  then  be  achieved  by  a 
shock  or  by  a  system  of  weaker  compression  waves.  The  boundary 
layer,  by  now  turbulent,  passes  through  recompression  without 
separating.  The  flow  accommodates  to  higher  a  by  gradual 
forward  movement  and  weakening  of  the  shock,  without  discontin- 
uities  in  Cn  or  Cm*  In  fully  developed  stall,  as  depicted  by 
the  data  for  a  =  16.6  degrees,  the  loading  has  reverted  to  the 
type  appearing  at  lower  Mach  numbers. 

The  main  features  of  the  static  loadincj  behavior  of  the  NACA 
0006  airfoil  are  closely  similar  to  those  of  the  Vertol  13006- 
.7,  but  the  sharper  leading  edge  and  the  lack  of  camber  lead 
to  their  development  at  lower  angles  of  attack. 

Dynamic  Effects  on  Loadings  During  Pitch  Oscillation 

Effect  of  Frequency 

Figure  20  presents  cycle  histories  of  Cn^  Cm^  and  five  selected 
chordwise  ACp's  for  the  Vertol  13006-.7  airfoil  oscillating  at 
12  Hertz  and  M  =  0.3.  The  data  shown  represent  the  average  of 
ten  consecutive  cycles  of  oscillation.  Steady  Cn  and  Cm  values 
corresponding  to  the  instantaneous  a  are  also  shown  for  compar¬ 
ison.  Note  that,  even  at  such  a  comparatively  low  frequency, 
a  significant  postponement  of  stall  is  evident  in  the  Cn  and 
Cm  traces. 
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Figure  21  shows  the  corresponding  sequence  of  chordwise  loading 
distributions  at  20-degree  intervals  throughout  the  cycle.  The 
first  loading,  at  an  a  before  stall  commences,  precisely  matches 
the  theoretical  steady-flow  loading.  This  is  true  of  all  cases 
where  stall  is  absent.  When  stall  does  occur,  as  at  o  «  20 
degrees,  it  is  initiated  by  collapse  of  the  leading-edge  ACp , 
and  this  effect  is  rapidly  transmitted  toward  the  trailing  edge. 
The  loadings  shown  between  0  =  80  and  G  =  140  degrees  are 
characteristic  of  fully  developed  stall.  Recovery  from  the 
stall  condition  begins  at  the  leading  edge,  with  a  return  of 
the  distribution  to  tho  theoretical  shape,  as  shown  for  0  =  140 
through  18J  degrees.  Tho  inception  of  stall  and  the  fully 
developed  stall  loadings  are  similar  for  all  Mach  numbers 
(0.2  to  0.6)  tested  at  this  frequency.  The  corresponding  data 
for  the  NACA  0006  behave  in  the  Seine  manner. 

The  effect  of  increasing  the  frequency  of  oscillation  to  68 
Hertz  is  shown  in  Figures  22  and  23,  which  arc  directly  compa¬ 
rable  to  Figures  20  and  21.  The  difference  in  dynamic  effects 
is  striking.  In  Figure  22  it  is  seen  that  a  marked  increase 
in  the  stall  delay  has  resulted  in  a  nearly  sinusoidal  Cm  trace, 
and  the  Cm  break  characteristic  of  the  lower  frequencies  no 
longer  exists.  Furthermore,  interesting  differences  in  the 
development  of  the  stalled  flow  are  indicated  by  the  loading 
data. 

The  sequence  of  chordwise  loading  plots  in  Figure  23  shows  that 
stall  begins  at  the  leading  edge.  A  wavelike  bulge  (presumably 
due  to  a  low-pressure  zone  on  the  upper  surface)  then  appears 
just  aft  of  the  leading  edge  and  moves  slowly  downstream.  (Its 
rate  is  about  0.1  chord  per  20-degree  interval.  At  this  re¬ 
duced  frequency,  the  free-stream  wind  velocity  corresponds  to 
0.44  chord  per  20-degree  interval.)  It  is  likely  that  the 
phenomenon  responsible  is  a  free  vortex  shed  from  the  leading 
edge  when  flow  separation  first  occurs,  and  which  then  drifts 
downstream  in  tlie  "wind  shadow"  of  the  airfoil. 

The  loadings  at  0  =  340  degrees  (i.e.,  -20  degrees)  ,  0  degrees, 
and  20  degrees  show  noticeable  deviations  from  the  steady-flow 
theoretical  line,  though  no  signs  of  stall  are  present.  These 
deviations  are  attributed  to  the  pitching  motion  which,  at 
this  point,  effectively  adds  camber  to  the  airfoil  by  producing 
a  relative  wind  component  in  the  downward  sense  ahead  of  the 
pitch  axis,  and  in  the  upward  sense  aft  of  it.  The  redistri¬ 
bution  of  loading  to  be  expected  therefrom  is  consistent  with 
the  noted  deviations. 

Effects  of  Mach  Number 


Figures  24,  25,  and  26  show  cycle  histories  of  ,  Cm,  and 
selected  ACp ' s  for  M  =  0.2,  0.4,  and  0.6,  respectively,  at  a 
nominal  reduced  frequency  of  0.2.  The  mean  angle  of  attack 
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Is  In  each  case  approximately  equal  to  that  for  steady-a  stall, 
while  Aa  -  5  degrees. 

The  force  and  moment  data  are  quite  similar  for  all  three  con¬ 
ditions.  In  each  case  goes  about  0.8  beyond  the  static 
level,  reaching  a  maximum  well  after  flow  separation  is  indi¬ 
cated  by  collapse  of  the  loading  at  the  leading  edge.  At  the 
same  time.  Cm  reaches  a  minimum  value  of  about  -0.3.  The 
recovery  process,  as  indicated  by  the  leading-edge  loading,  is 
compT  3te  at  the  220-degree  point  in  the  cycle.  Between  200 
degrees  and  260  degrees,  it  appears  that  dynamic  loading  due 
to  nose-down  pitch  rate  reduces  Cn  and  produces  a  positive  Cm. 

Substantial  differences  are  visible  only  in  the  loading  data, 
mainly  at  the  beginning  and  end  of  the  cycle.  The  meucimum 
local  loading  is  strongly  limited  by  Mach  number,  declining 
from  ACp  ■  7.0  atM“  0.2  to  ACp  =  4.5  at  M  =  0.6.  Transonic 
flow  effects  account  for  this  limitation  effect,  as  noted  pre¬ 
viously. 

CYCLE-TO-CYCLE  FORCE  AND  MOMENT  VARIATIONS 

All  data  presented  elsewhere  in  this  report  were  obtained  by 
averaging  the  measurements  taken  from  groups  of  ten  or  more 
consecutive  cycles  of  oscillation.  The  dynamic  response  of  a 
rotor  blade  will,  of  course,  depend  on  the  loads  acting  over 
individual  cycles.  The  question  of  variation  from  cycle  to 
cycle  is  therefore  of  practical  importance. 

Figure  27  shows  typical  Cn  and  Cm  traces  obtained  during  four 
successive  cycles  of  pitch  oscillation  for  the  Vertol  13006-.7 
airfoil  where  strong  stall  effects  are  present.  The  variations 
in  both  Cn  and  Cm  negligible  where  stall  effects  are  absent, 
as  shown  by  that  part  of  the  cycle  between  stall  recovery  and 
inception.  During  the  initial  period  of  stall,  Cn  and  espe¬ 
cially  Cm  show  very  repeatable  behavior.  During  the  stall 
recovery  period,  however,  a  wide  variety  of  patterns  may  be 
followed. 

Since  Cn  and  Cm  behave  in  a  regular  manner  when  attached  flow 
prevails  and  also  when  the  flow  breakdown  begins,  it  may  be 
inferred  that  the  marked  scatter  in  the  stalled-flow  behavior 
is  not  caused  by  wind  tunnel  flow  turbulence  or  by  variations 
in  the  model  airfoil  motion.  Instead,  these  flow  variations 
during  stall  recovery  are  attributed  to  the  presence  of  large- 
scale  turbulent  eddies  formed  above  the  airfoil  as  suggested 
from  the  previous  discussion  of  dynamic  loadings  during  stall. 
The  random  nature  of  turbulent  flow  is  well  known  and  can  be 
expected  to  result  in  loading  variations  susceptible  only  to 
statistical  description. 
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DRAG  IN  PITCHING  OSCILLATION 


Because  of  difficulties  In  measuring  the  instantaneous  airfoil 
drag  under  unsteady-flow  conditions,  as  reported  in  Reference 
2,  only  time-average  drag  values  were  determined.  The  average 
drag  in  pitch  oscillation  was  obtained  by  using  a  slowly  tra¬ 
versing  pitot-static  probe  to  survey  the  wake  momentum  profile. 

Drag  data  were  obtained  for  both  airfoils.  Measurements  for 
the  NACA  0006  airfoil,  teOcen  at  M  «  o.4,  are  compared  in  Figure 
28  to  the  steady-flow  Cd,  which  is  shown  as  a  dashed  line.  The 
oscillatory  drag  data  agree  with  the  trends  for  thicker  airfoils 
given  by  Reference  2  and  show  that  the  average  Cj)  in  pitch 
oscillation  is  always  higher  than  the  steady  level,  increasing 
with  frequency  between  12  and  48  Hertz.  A  small  part  of  this 
increase  may  be  due  to  the  fact  that  Aa  increases  with  frequen¬ 
cy  because  of  mechanical  effects. 

Figure  29  shows  drag  data  obtained  for  the  Vertol  13006-.7 
airfoil  at  M  -  0.4  and  0.6.  At  0.4,  the  data  trends  are 
similar  to  those  for  the  symmetrical  airfoil  but  show  a  2-degree 
advantage  in  the  angle  of  attack  for  the  drag  increase  due  to 
stall,  which  is  consistent  with  the  static  behavior.  At  M  > 

0.6,  oscillatory  data  were  obtained  only  at  a  frequer.cy  of  12 
Hertz. 

COMPARISON  WITH  PREVIOUS  RESULTS 


Static  and  dyneunic  data  from  this  test  on  the  NACA  0006  and 
Vertol  13006-.7  airfoils  are  compared  in  Figures  30  through 
33  with  the  results  of  a  similar  series  of  tests  on  the  symmet¬ 
rical  NACA  0012  (modified)  and  the  cambered  Vertol  23010-1.58 
sections. 

The  steady-flow  Cn  and  Cm  characteristics  are  compared  at  M  » 
0.4  in  Figure  30.  Both  thick  airfoils  show  the  large,  eUorupt 
changes  in  Cn  and  Cm  characteristic  of  leading-edge  stall^. 

This  stall  is  caused  by  a  sudden  separation  of  the  flow  over 
the  entire  upper  surface  of  the  airfoil.  No  appreciable 
buffeting  was  observed  from  the  pressure  traces  for  these  air¬ 
foils  below  stall.  The  thin  airfoils,  on  the  other  hand,  show 
a  change  in  the  lift-curve  slope  but  little  or  no  loss  of  lift, 
and  a  comparatively  gradual  development  of  the  nose-down 
pitching  moment  due  to  stall.  (This  is  characteristic  of 
thin-airfoil  stall,  as  noted  above.)  Buffeting,  detected  from 
ACp  fluctuations  on  the  analog  records,  begins  before  the  lift 
curve  changes  slope. 

Comparisons  of  the  characteristics  of  the  dynamic  data  can  be 
summarized  in  terms  of  two  important  parameters:  the  dynamic 
Cnmax  frequencies  corresponding  to  a  typical  helicopter 
rotor  1/rev,  and  cycle  damping  at  the  rotor  blade  first 
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torsional  natural  frequency.  For  convenience,  the  cycle!  damp¬ 
ing  is  divided  by  the  theoretical  damping,  based  on  Theodorson's 
formulation^ . 

Curves  of  Cnm;^X  ^'g^^sus  Mach  number  are  presented  in  Figure  31. 
Cycle  damping  of  the  four  airfoils  is  compared  in  Figure  32  for 
6/rev  pitch  motion  (k  =  0.74)  at  M  =  0.2.  The  cambered  air¬ 
foils  show  a  substantial  advantage  over  the  symmetrical  air¬ 
foils,  both  in  the  inception  and  in  the  magnitude  of  negative 
damping  at  a  given  u. 

The  region  of  instability  (based  on  zero  damping)  is  shown  in 
Figure  33  for  an  equivalent  frequency  of  6/rev.  Both  camber 
and  airfoil  thickness  are  seen  to  postpone  the  regions  of 
negative  damping  to  higher  angles  of  attack. 

The  general  reduction  in  the  range  of  instability  with  in¬ 
creasing  Mach  number  is  due  to  transonic  flow  effects.  For  the 
thicker  airfoils,  stall  caused  by  the  sudden  flow  separation 
on  the  upper  surface  at  Mach  numbers  up  to  0.5  is  replaced  by 
a  stable  expansion-fan  and  shock  system  at  higher  M.  Separa¬ 
tion,  if  present,  occurs  aft  of  the  shock. 

Transonic  flow  effects  appear  sooner  in  the  case  of  6-percent- 
thickness  airfoils,  but  they  are  limited  to  smaller  regions 
close  to  the  leading  edge  because  of  the  sharper  nose  profiles. 

COMPARI./JN  WITH  THEORETICAL  PREDICTIONS 


In  Figure  34,  chordwise  loading  distributions  under  steady-flow 
conditions  below  stall  are  compared  with  theoretical  predictions 
for  both  airfoils  at  M  =  0.2,  where  compressibility  effects  are 
virtually  absent.  Data  for  the  NACA  0006  were  taken  directly 
from  Reference  6,  while  those  for  the  Vertol  13006-.7  airfoil 
were  developed  using  the  vorticity  polygon  method  of  Reference 
7.  Agreement  is  excellent  in  both  cases. 

The  experimental  dynamic  force  and  moment  derivatives  for  both 
airfoils  at  small  values  (no  stall  present)  are  compared 
with  theoretical] y  predicted  data  in  Figure  35.  The  theoreti¬ 
cal  values,  based  on  thin  airfoil  representation  with  the 
effects  of  compressibility  and  tunnel  walls  (method  of  images) 
included,  were  calculated  with  a  computer  program  supplied  by 
NASA®.  The  data  show  essentially  identical  characteristics 
for  each  airfoil.  The  differences  between  test  and  theory  are 
similar  to  the  trends  noted  for  thicker  airfoils,  as  reported 
in  Reference  2. 

The  experimental  values  of  normal  force  amplitude  are  slightly 
lower  than  predicted.  This  is  consistent  with  the  universally 
noted  reduction  in  static  lift-curve  slope  caused  by  viscosity. 
The  trend  to  higher  lift-curve  slope  as  Mach  number  is  increased 
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is  in  good  agreement  with  the  predictions.  The  data  for  lift 
phase  agree  very  well  with  the  theory,  maximum  differences 
being  less  than  5  degrees. 

The  pitching  moment  amplitudes  agree  well  in  trends  with  Mach 
number  but  show  a  reduction  in  magnitude,  similar  to  that 
noted  for  C^.  This  reduction  is  attributed  to  viscous  effects. 
The  Cm  phase  data  show  appreciable  departure  from  predictions 
away  from  the  midranges  of  reduced  frequency,  with  maximum 
differences  up  to  30  degrees.  These  variations  result  from 
the  extreme  sensitivity  of  the  pitching  moment  phase  to  the 
pitch  reference  center,  coupled  with  light  shifts  of  the 
aerodynamic  center  from  the  quarter-chord  point  for  the  addi¬ 
tional  loading  due  to  oscillation. 
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CONCLUSIONS 


1.  Thin  airfoils  oscillating  in  pitch,  at  the  eunplitudes 
and  frequencies  typically  experienced  by  blade  elements 
of  helicopter  rotors,  show  a  strong  dynamic  postponement 
of  lift  and  moment  stall  to  higher  angles  of  attack  than 
experienced  under  steady-flow  conditions. 

2.  Substantial  increases  in  the  maximum  normal  force  edjove 
steady-flow  values  are  attained  durj.ng  oscillation,  and 
Increase  with  frequency  but  decline  as  the  Mach  number 
increases  to  0.6.  Increases  in  the  meucimum  normal  force 
capability  are  about  equal  to  those  attained  by  the 
thicker  airfoil  sections  reported  in  Reference  2. 

3.  The  Vertol  13006-.7  airfoil  stalls  at  higher  amgles  of 
attack  than  the  NACA  0006  section.  This  is  consistent 
with  the  differences  in  steady-flow  behavior  and  is  un¬ 
affected  by  oscillation  frequency. 

4.  Pitch  oscillations  at  angles  of  attack  where  stall  takes 
place  can  result  in  net  energy  inputs  from  the  airflow 
into  the  airfoil  system.  These  conditions  of  negative 
damping  imply  the  possibility  of  divergent  torsional 
oscillations  in  a  rotor  blade.  The  inception  of  nega¬ 
tive  damping  is  shifted  to  higher  Oo  by  increasing 
frequency,  but  it  occurs  earlier  as  Mach  number  is  in¬ 
creased.  At  frequencies  between  4/rev  and  6/rev,  cor¬ 
responding  to  a  typical  rotor  blade  first  torsional 
vibration  mode,  the  width  and  Intensity  of  negatively 
damped  regions  diminish  at  the  higher  Mach  numbers.  At 
M  =  0.6,  negative  damping  at  4/rev  frequency  has  dis¬ 
appeared  but  .8  present  at  6/rev. 

5.  The  flow  processes  governing  stall  development  and  recov¬ 
ery  are  initiated  by  conditions  at  the  leading  edge  and 
rapidly  progress  toward  the  trailing  edge.  The  character 
of  the  chordwlse  loading  during  stall  and  reattachment 
shows  considerable  variation  over  the  range  of  test  con¬ 
ditions  studied.  It  indicates  that  large-scale  vortex 
shedding  is  Induced  on  the  airfoil  upper  surface,  and 
that  its  intensity  is  strongly  dependent  on  the  reduced 
frequency. 
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RECOMMENDATIONS 

1.  The  data  presented  In  this  report  should  be  further  ana¬ 
lyzed  with  a  view  to  the  establishment  of  a  mathematical 
description  of  dyn2unic  stall.  An  empirical  formulation 
of  this  type  is  a  necessary  first  step  in  the  practical 
application  of  these  data  to  rotor  dynamic  and  aerody¬ 
namic  computer  analyses. 

2.  A  program  of  analysis  should  be  undertaken  to  develop  a 
detailed  understemdlng  of  the  flow  processes  affecting 
dynamic  stall,  such  as  local  supersonic  flow  emd  vortex 
shedding.  The  program  would  require  supplemental  experi¬ 
mental  work  in  flow  visualization  and  the  measurement  of 
absolute  pressures  to  define  local-flow  Mach  numbers  on 
the  airfoil  surfaces. 
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Figure  5.  Effect  of  Mean  Angle  of  Attack  on  Dynamic  Cjj 
and  Cm  Versus  o  for  Pitch  Oscillation  of  the 
Vertol  13006-.7  Airfoil  at  M  =  0.2,  f  =  12  Hertz, 
and  Aa  =  5*. 


25 


LEGEND 


LINE 

TYPE 

Aa 

f 

k 

0* 

TP 

9.»J 

12.34 

0.108 

10.07 

8107.3 

7.01 

11.77 

0.103 

12.48 

8097.4 

4.S4 

11.92 

0.115 

14.82 

8018.3 

PITCHING  MOMENT  COEFFICIENT  (Cm)  NORMAL  FORCE  COEFFICIENT  (Cm) 


ANGLE  OF  ATTACK  (a)  -  DEGREES 


ANGLE  OF  ATTAC 


Figure  8.  Effect  of  Mach  Number  on  the  Dynamic  Cfj  and 

for  the  Vertol  13006-.7  Airfoil  at  Aa  =  5°.  ” 


31 


ANGLE  or  ATTACK  (o)  -  DEGREES 


ANGLE  OF  ATTACK  (a)  -  DEGREES 


namic  Cm  and  C. 
at  Ao  =  5».  ^ 


PITOIIMG  NamT  COBPFICIIHT  <0^)  NOHUL  POMCB  COBmCZINT  (0^) 


TT  f  (Pi)  k  N  ia  a, 

S034.4  49.17  0.331  0.404  S.58  9.M 


0  9  10  19  20  29 

or  Maa  (•)  -  DMHU 


Figure  10.  Effect  of  Leading-Edge  Design  on  the  Dyn2unic 
Cn  and  at  f  =  48  Hertz,  M  =  0.4,  and 
Aa  s*  5® . 
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MACH  NUMBER  (M) 

Figure  11.  Summary  of  Dynamic  Maxima  and  Minima  for 
Pitch  Oscillation  at  Aa  *  5®  and  f  =  12  Hertz. 
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Figure  12.  Marimum  Cn  Attained  for  the  Vertol  13006-.7 
Airfoil  in  Pitch  Oscillation  at  Aa  =  5°. 
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Figure  13.  Maximum  Cm  Attained  for  the  NACA  0006 

Airfoil  in  Pitch  Oscillation  at  A.'i  =  5°. 
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DAMPING/THEORETICAL  DAMPING 


Figure  14,  Effect  of  Oscillation  Frequency  on  Pitch 
Damping  Characteristics  for  the  Vertol 
13006-.7  Airfoil  at  M  =  0.2  and  Aa  =  5®. 
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Figure  16.  Effect  of  Oscillatory  Amplitude  on  Pitch  Damping 
for  the  Vertol  13006-.7  Airfoil  at  M  =  0.4  and 
f  =  48  Hertz. 
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Figure  18.  Comparison  of  Pitch  Damping  Stability  Limits 
for  the  Vertol  13006-.7  and  NACA  0006 
Airfoils  at  4/Rev  Frequency. 
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PITCHING  NOH£NT  COEFFICIENT  (C>,)  NORMAL  FORCE  COEFFICIEtlT  (C||) 


2.S 


AIRFOIL  PITCH  REFERENCE  (6)  -  DECREES 


Figure  20.  Cycle  History  of  C^,  Cm,  and  ACp  for  the  Vertol 
13006-.7  Airfoil  Oscillating  in  Pitch  at  M  =  0.3 
f  =  12  Hertz,  Aa  =  5®,  and  =  10°.  ' 
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Cycle  History  of  Cn#  Cm*  and  ACp  for  the  Vertol 
13006-.7  Airfoil  Oscillating  in  Pitch  at  M  >  0. 
f  ■  72  Hertz /  Aa  »  5®,  and  a*  ■  10®. 
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Pigur«  27.  Successive  Cycle-to-Cycle  Variations  in  Cn  and 
Cm  During  Pitch  Oscillation  for  the  Vertol 
13006-.7  Airfoil. 
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Figure  30.  Comparison  of  Steady-Flo«r  Cn  and  Cm  Characteristics 
With  Data  of  Reference  2  at  M  =  0.4. 
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BEST  DYNAMIC 


Figure  31.  Comparison  of  Maximum  Attainable  Normal  Force 
Characteristics  With  Data  of  Reference  2,  at 
f  =  1/Rev  Frequency  and  La  =  5®. 
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Figure  32.  Comparison  of  Aerodynamic  Damping  With  Data  of 
Reference  2,  at  6/Rev  Frequency,  M  =  0.4,  and 
Aa  =  5“. 
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Figure  33.  Comparison  of  Aerodynamic  Instability  Regions 
With  Data  of  Reference  2,  at  6/Rev  Frequency 
and  A a  =  5®. 
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APPENDIX  I 

STATIC  (NONOSCILLATORY)  TESTS 


Each  airfoil  was  tested  under  the  steady-a  condition  over  the 
full  range  of  a  and  M  at  which  dyneimic  data  were  obtained. 

Figures  36  through  38  show  steady  and  versus  a  charac¬ 
teristics  for  both  airfoils.  The  position  of  sharp  changes 
in  the  airfoil  characteristics  was  verified  from  continuous 
emalog  records . 

The  cambered  airfoil  has  higher  Cni4;vx  o>STALL  than  the 
symmetrical  section.  Neither  profile  is  subject  to  large 
discontinuities  on  the  Cn  curve.  Below  M  =  0.4,  both  airfoils 
exhibit  a  leading-edge  laminar  separation  bubble  which  is  a 
well-known  characteristic  of  thin  airfoils^.  The  onset  of 
this  phencxnenon  is  indicated  by  the  notch  in  the  Cn  curve, 
which  occurs  later  with  the  cambered  airfoil  and  precedes  full 
stall  development. 

Both  airfoils  show  well-behaved  pitching  moment  behavior  at 
angles  of  attack  below  stall.  The  presence  of  small  changes 
in  dCM/da  with  increasing  Mach  number  for  the  NACA  0006  air¬ 
foil  reflects  a  slight  forward  shift  of  the  aerodyn£unic  center. 
Stalling  characteristics  indicated  consistent  trends  for  the 
range  of  Mach  numbers  tested,  without  the  abrupt  drop  shown 
by  the  thicker  airfoils  in  Reference  2. 

In  Figures  39  through  41,  Cn  versus  Cm  curves  are  presented 
for  each  airfoil.  Comparison  with  unpublished  data  obtained 
during  1965  from  tests  using  the  same  wind  tunnel  facilities 
is  given  for  M  «  o.4  emd  0.6  in  Figure  42.  Good  agreement 
was  obtained  for  both  airfoils. 

Static  drag  data  from  wake  traverse  measurements  are  shown  in 
Figure  43  for  both  airfoils.  The  tests  on  the  NACA  0006, 
limited  to  N  »  0.4,  show  a  lower  a  for  drag  rise  due  to  stall. 
This  is  consistent  with  the  normal  force  and  moment  data. 

Figure  44  presents  static  Cn  versus  Cq  data  for  the  Vertol 
13006-.7  airfoil  in  comparison  with  the  previous  measurements. 
Good  agreement  is  obtained  at  both  M  ^  o.4  and  0.6. 
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Figure  36.  Static  C[|  and  Cm  Vlarsos  a  Characteristics  for 
the  Vertol  13006-.7  Airfoil  at  M  =  0.2.  0.4. 
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PITCHING  NONENT  COEFFICIENT  (C^^) 

Figure  39.  Static  Cn  Versus  CH  Characteristics  for  the 

Vertol  13006-.7  Airfoil  at  M  -  0.2,  0.4,  and  0.6. 
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PITCHING  MOMENT  COEFFICIENT  (C„) 

M 


Static  Cn  Versus  Cj^  Characteristics  for  the  NACA 
0006  Airfoil  at  M  =  0.2,  0.3,  0.4,  and  0.6. 
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APPENDIX  II 

WIND  TUNNEL  WALL  CORRECTIONS 

Wind  tunnel  boundary  corrections  were  developed  by  using  a  com¬ 
puter  program  provided  by  NASA  (based  on  the  formulation  in 
References  9  and  10)  to  calculate  CWf  Cm^  and  phase  for  an  air¬ 
foil  oscillating  in  a  wind  tunnel.  Comparative  free-air  values 
were  obtained  from  the  data  teibles  in  Reference  11.  The  theo¬ 
retical  values  were  based  on  compressible  thin  airfoil  theory, 
with  the  method  of  images  used  to  simulate  the  tunnel  walls. 

The  amplitude  ratios  of  Cn  and  Cm  in  free  air  to  the  values  in 
a  tunnel,  as  well  as  phase  differences,  are  shown  in  Figure  45. 
The  corrections  to  magnitude,  which  are  largest  at  low  values 
of  k  and  high  Mach  number,  are  all  less  than  14  percent.  The 
corrections  to  phase  are  all  less  than  10  degrees,  except  for 
the  pitching  moment  at  M  =  0.6  and  low  values  of  k. 

The  corrections  calculated  here  were  found  to  be  small  in  com¬ 
parison  with  the  large  changes  in  both  magnitude  and  phase  due 
to  stall,  and  therefore  they  have  not  been  included  in  the 
wind  tunnel  results. 
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VJall  Corrections  for  Pitching 
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111.  lOONlOIIIM*  UILITARV  ACTIVITY 


Volume  I  of  a  2-volumj  report 


U.S.  Army  Aviation  Materiel 
Laboratories 
Fort  Eustis,  Virginia 


AoiTRACv  report  presents  the  results  of  an  experimental  investiga- 

tion  of  rotor  blade  dynamic  stall.  Forces  and  moments  in  two* 
dimensional  flow  were  determined  for  two  thin  airfoil  sections  (NACA 
0006  and  Vertol  13006*. 7)  usually  considered  in  advanced  helicopter 
rotor  concepts  by  measuring  differential  pressure  during  oscillatory 
pitch  motions.  The  Mach  number,  Reynolds  number,  and  angle  of  attack 
prevailing  in  the  retreating-blade  stall  region  were  investigated  up 
to  typical  first* torsion  mode  natural  frequencies. 

Pitch  oscillation  was  shown  to  increase  the  angle  of  attack  at  which 
stall  occurred,  with  a  corresponding  increase  in  maximum  normal  force. 
Hysteresis  effects  on  the  pitching  moment  response  produced  regions 
of  negative  aerodynamic  damping  in  oscillations  at  mean  angles  of 
attack  near  the  steady-flow  stall  values  for  all  Mach  numbers  from 
0.2  to  0.6. 

It  is  shown  that  dynamic  stall  effects  are  strongly  related  to  the 
reduced  frequency  parameter.  Significant  evidence  of  compressibility 
effects  appears  only  at  Mach  numbers  higher  than  0.4.  At  lower  Mach 
numbers,  the  data  indicate  that  the  two  effects  of  decreased  Mach 
number  and  airfoil  camber  are  fully  comparable  in  terms  of  dynamic 
stall  patterns  and  increase  both  the  usable  angle  of  attack  and  the 
maximum  normal  force. 


